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Summary. — At high nuclear density and small temperature due to the asymptotic 
freedom property of Quantum ChromoDynamics and to the existence of an attractive 
channel in the color interaction, diquark condensates might be formed. Since these 
condensates break the color gauge symmetry, this phenomenon has the name of 
color superconductivity. In the last few years this has become a very active field of 
research. While a direct experimental test is still missing, color superconductivity 
might have implications in astrophysics because for some compact stars, e.g. pulsars, 
the baryon densities necessary for color superconductivity can probably be reached. 



1. Nuclear matter and QCD 

At small nuclear density and low temperature the relevant degrees of freedom of 
nuclear matter are the hadrons: nucleons, hyperons, mesons. At high density and/or 
high temperature the fundamental fields are those of quarks, spin 1/2 particles carrying 
electric charge and a quantum number called color, and those of gluons, massless spin 
1 particles carrying color in eight possible different ways. The corresponding theory 
is called Quantum ChromoDynamics (QCD). It resembles Quantum Electrodynamics 
(QED), where the role of the fundamental fermion is played by the electron and that of 
the eight gluons by the photon (carrying no charge and therefore unique). Similarly to 
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QED, QCD is a gauge theory based on a local unitary gauge group. For QCD it is SU(3), 
while QED is based on the abelian group U(l). Quarks belong to the fundamental and 
gluons to the adjoint representation of SU(3), which is why each quark has three and 
each gluon eight possible colors. Besides color, quarks carry another quantum number, 
flavor, that affects their mass. Six different flavors are known, but we will consider here 
only the lightest quarks, called up, down and strange (with m u ~ nid <C m s ). A peculiar 
property of QCD is Asymptotic Freedom (AF), which means that at large momenta the 
effective QCD coupling constant becomes smaller and smaller with increasing momenta. 

In the ordinary world, i.e. at zero temperature and zero density, the study of the QCD 
degrees of freedom is made difficult by the fact that quarks and gluons are confined inside 
hadrons. Studying the QCD vacuum state under extreme conditions can therefore help 
our understanding of strong interactions, because their description is simpler in these 
regimes, due to the AF property of QCD. In particular one might perform the study in 
two very different regimes: high temperature and small density or small temperature and 
high density. The first regime is presently studied at the Relativistic Heavy Ion Collider 
at Brookhaven, USA, and will be investigated in the near future at the Large Hadron 
Collider at CERN. The other regime is not accessible for the present at the accelerators, 
but can be indirectly investigated by the study of the properties of compact stars such 
as pulsars, where both conditions, small temperatures and high densities, are likely to 
occur. 

As we noted quarks are spin 1/2 fermions with three colors, conventionally called 
red, green, blue (r,g,b) and, if light, in three flavors (u,d,s). Since they belong to the 
representation 3 of the color symmetry group, when two quarks scatter they are either in a 
sextet (color symmetric) or in antitriplet (color antisymmetric) state. QCD interactions 
are mediated by one gluon exchange (in the perturbative regime). In the antitriplet 
channel this interaction is attractive. As a consequence of the Cooper's theorem the 
perturbative vacuum is unstable and one expects the formation of Cooper pairs of quarks 
(diquarks) and color superconductivity (this is so because the condensate is colored, not 
white as in the quark-antiquark channel, relevant at zero density). Moreover the two 
quarks scatter preferably in the spin channel, since this corresponds to a larger phase 
space. At extreme densities, considering QCD with 3 flavors, the attractive interaction 
in color antisymmetric channel, together with 5 = implies, by the Pauli principle, 
antisymmetry in flavor. This is the basic idea of the models that we discuss in the next 
section. 

2. The true vacuum 

Since the naive vacuum is unstable we should try to find the true superconducting 
vacuum state. To do that we have to show that the formation of a quark-quark con- 
densate occurs because of the color interaction. The interaction term is fundamentally 
provided by gluon exchanges and to compute it and derive the gap equation one can 
follow different approaches, such as, for example, the one of Bailin and Love in based 
on the resummation of bubble diagrams (for previous work see references in £Q). For this 
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method and similar approaches we refer to the reviews 0. 

For pedagogical purposes it is however sufficient to consider a Nambu-Jona Lasinio 
[3] (NJL) interaction (4 fermion interaction) that mimics the QCD interaction. In this 
approximation the hamiltonian is: 

(1) H = j 'd 3 x($(i$ + li 1o )tP + G^T A ^^T a ^) . 

Here G is a coupling constant with dimensions (mass) -2 and T A are SU(3) generators. 
Notice that a form factor is included in the definition, to take into account the asymptotic 
freedom property of QCD. It would act as an ultraviolet cutoff. 

In order to obtain the gap equation we make a mean field ansatz: 

(2) - (t'Cfe). 

C is the charge conjugation matrix. Here and below Greek indices a, [3 are color indices, 
Latin indices i, j are flavor indices. Let us now show that one expects for the vacuum 
expectation value of the quark bilinear the result 

(3) {i>ai C ^Pi) = tafiieijK^T- 

(sum over repeated indices). It is useful to note that, by evaluating in mean field 
approximation, i.e. substituting to a pair of fermion fields in their mean value given 
by (J3J and J3J), the interaction term takes the form 

(4) ~ A^V + h.c. , 

i.e. the fermions acquire a Majorana mass. 

To prove ® we note that, as we have already discussed, QCD interaction favors 
antisymmetry in color; the pair has zero angular momentum (in this way the entire 
Fermi surface is available and the effect becomes macroscopic); if it has also spin zero 
the fermion pair must be in an antisymmetric state of flavor to produce an antisymmetric 
wavefunction. Neglecting the sextet contribution, which is small, leaves us with J3J). 

Clearly the form of the matrix A/#- corresponding to the true vacuum state depends 
on dynamical effects. In the literature various cases have been discussed: 

1. 2 flavor Superconductivity (2SC): A IK = A5i 3 S K3 g] ; 

2. Color-Flavor-Locking (CFL): A IK = S IK A ; 

3. gapless Color-Flavor-Locking (gCFL): Aik — SikAk 013 ; 

4. 2 flavor LOFF phase: A IK = 5 n 5 K ?,A{v) El EB CI] ! 

5. 3 flavor LOFF phase:A /isr = S IK A K (r) \T2\ . 
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It is an incomplete list, for other results see 23,23] and references in 

In the case of the LOFF phases the r dependence of the condensates produce an in- 
homogeneous type of superconductivity. Various r dependences have been discussed, e.g. 
the single plane wave or the sum of several plane waves (see below) . In the cases of the 
gCFL phase and the 3 flavor LOFF phase one has to take into account the differences 
of quark masses (as we noted above, m s is significantly larger that the up and down 
masses) and the differences of quark chemical potentials (jj, u , /x^, /i s are in general differ- 
ent). This introduces more parameters in the description, but their number is reduced 
by imposing that the quark matter is electrically and color neutral (such constraint is 
automatically implemented in the CFL model). We will concentrate our attention on the 
above mentioned models. 

3. Homogeneous color superconductivity 

To show that color superconductivity can exist, let us consider the first two cases (2SC 
and CFL). The 2SC model corresponds to the presence in the condensate of only the u 
and d flavors: this means that this ansatz may be relevant in the situation characterized 
by an intermediate chemical potential m Ul <C [i <C m s . Therefore one can imagine 
of having taken the limit m s — > oo and the theory is effectively a two-flavor theory. 
The CFL model holds in the other limit m u , nid, m s <C M- It takes its name from the 
property 

3 

(5) ^ £ a f3K£ijK = S a iSpj — 6 a jSf3i , 

K=l 

by which the color and flavor indices are locked. 

Let us now introduce a,(V ,b, ... annihilation/creation operators of particles and holes: 

(6) i> ai = Y,u k (a k e ikx + bte- ik3: ) . 

k 

The existence of the condensate gives rise to the following hamiltonian 
H = H i} k ~ ^ a t ak + ( fc + V)b{b k ) 

k 

(7) + ]T j e-^aua-u + je+ l %{bl k + he 

fc 

where we have omitted for simplicity any color-flavor indices as we wish to stress the 
mechanism that produces the condensate. Let us now perform a unitary Bogoliubov 
transformation. To do that we introduce annihilation and creation operators for quasi- 
particle and quasiholes 

j/fc = cos Oak — e J * sin^a^ , 
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(8) zu — cos tpbk — e l * sin <pbl_ 



fc • 



By an appropriate choice of the parameters of the transformation we can transform the 
original hamiltonian in a new hamiltonian describing a gas of non interacting quasipar- 
ticles 

(9) H = ^2[e v (k)y{y k + e z (k)4z k 

with 

e y (k) = y/(k - m) 2 + A 2 , 



(10) e z {k) = y/(k + /i) 2 + A 2 . 

These equations show two effects: 

1. The quasiparticles and quasiholes are free; 

2. in the dispersion laws of the quasiparticles a mass term proportional to A (the gap 
parameter) appears. 

To obtain (J3J, i.e. a free hamiltonian, the parameters have to be chosen as follows: 

cos2fl= - , 

(11) cos2^- fc + A< 



vRfc + m) 2 + a 2 

We note that, differently from the original a*,, bk operators that annihilate the false 
vacuum, the quasiparticles annihilation operators destroy the true vacuum 

(12) 2/fc|0>=0, z fc |0>=0. 

We have still to prove that A ^ 0. We do that and get an equation for A by 
substituting for yk,Zk in cq. (J3J. We get in this way an integral equation that has the 
following schematic form (gap equation): 



(13) A = C / d 3 k 



A 



V(A - M) 2 + A 2 v/(fc + /x) 2 + A 2 

We see immediately that the origin of the instability of the false vacuum lies in the first of 
the two terms in the r.h.s. of l|13|l : If A = there is no compensation for the divergence 
at k = fj,. 

To be more quantitative let us consider the two models in more detail. In the CFL 
model (m s = m u = rrid = 0) all the 3x3 = 9 quarks acquire a Majorana mass. The CFL 



6 



G. Nardulli 



condition gives two different set of eigenvalues. The first one comprises 8 degenerate 
masses Ai = A2 = ... = Ag and the second set the non degenerate mass Ao. The actual 
values of the gaps depends on the model and the approximations involved; typical values 
are 



for n « 400 - 500 MeV. 

Let us observe that, since we start with a massless theory, the QCD hamiltonian has, 
besides the local color symmetry SU(3), also a global chiral symmetry ST/ (3) l x SU(3)r. 
The pairing occurs for left-handed and right handed quarks separately. Therefore these 
condensates break SU(3) c ®SU(3)l®SU(3)r; finally, since the pair has baryonic number 
2/3 (each quark has baryonic number equal to 1/3 since the nucleon is made up by 
three quarks), also U(1)b is broken. It should be observed, however, that a diagonal 
SU(3) c +l+r symmetry remains unbroken. 

In the 2SC model (m s >/j> m u , m^) there are 4 massless quarks (the up and down 
quarks with colors 1 and 2), while the strange quark with any color and the u and d 
quarks with color 3 remain ungapped. Numerical results for this case are qualitatively 
in agreement with the results ((HI) of the CFL model. 

It is interesting to observe that the spontaneous breaking of global symmetries implies 
the existence of Nambu Goldstone bosons; for internal symmetries there are as many 
NGBs as there are broken generators G a . As they are massless, they are the lowest 
energy quasiparticles of the effective theory. On the other hand, for broken local gauge 
symmetries, by the Higgs- Anderson mechanism the gauge bosons acquire masses and 
there are no NGBs. In the CFL model the diquark condensate breaks SU (3) l X SU (3)r x 
SU(3) C x U(1)b to SU(3) c+ l+r x Z 2 . All the 9 quarks are massive and they belong to 
a SU(3) c+ l+r, singlet and a SU(3) c+ l+r octet. All the 8 gluons are massive and are 
degenerate. There are 8+1 Nambu-Goldstone bosons. In the 2SC model the condensate 
breaks SU(3) C <g> SU(2) L (g> SU(2) R X U(1) B down to SU(2) C x SU(2) L xSU{2) R . While 
the chiral group is unbroken, SU(3) C is broken to SU(2) C . Therefore 3 gluons remain 
massless and 5 acquire a mass. As to the other quasiparticles, there is one would-be NGB 
associated to the breaking of the axial symmetry; moreover of the 6 quarks (2 flavors in 
3 colors) 4 are massive and 2 are massless. 

One can show that, for arbitrarily large /i, the CFL model is favored ^HJ; however, 
for intermediate /1 other phases should be considered. At lower densities, in a range 
presumably more relevant for the study of compact stars, taking into account strange 
quark mass effects implies that a new phase, called gCFL (gapless CFL) phase, is favored. 
Let us therefore consider now this new phase. 

The Lagrangian for gluons and ungapped quarks with m u = rrid = and m s ^ can 
be written as follows (color, flavor and spin indices suppressed): 



(14) 



Ax 



A 8 = A w 20 - 100 MeV 



A « -2A , 



(15) 



£ = tp M + ^70) V 
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where M = diag(0, 0, m s ) is the mass matrix in flavor space and the matrix of chemical 
potential is given by 

(16) atf = (^6% - mQa) s aP + ^ 3 Tf + -J^if) 

= 1,3 flavor indices; a, (3 = 1,3 color indices). Moreover T 3 = |diag(l, — 1, 0), 
T 8 = j^diag(l, 1, —2) in color space and Q = diag(2/3, —1/3, —1/3) in flavor space; 
is the electrostatic chemical potential; /i3, fig are the color chemical potentials associated 
respectively to the color charges T3 and Tg; \Xh is quark chemical potential. 

Working in the mean field approximation and neglecting the antiquark (hole) con- 
tribution, one can derive the quasiparticle dispersion law, whose knowledge allows the 
evaluation of the grand potential. At zero temperature it is given by 

(17) n = tdpp 2 5> j (p)| + l^ + Af + A 2 .) 



i=i 

where A is the ultraviolet cutoff, ej (p) are the quasi particle dispersion laws and G is the 
Nambu-Jona Lasinio coupling constant. 

In any realistic astrophysical application of these models one has to impose electric 
neutrality (the global electric charge must vanish) and color neutrality (the color charges 
must vanish as well). In absence of these conditions, electric and color forces would 
destroy quark matter stability. In order to enforce electrical and color neutrality one has 
to minimize the grand potential with respect to hq , ^ and /is ■ Including the stationary 
conditions with respect to the gap parameters Ai,A2,A3 one ends up with a system 
of six equations which must be solved simultaneously. Once this system of equations 
is solved one may express Ai,A2,A3 and the chemical potentials fiQ,/J-3 and /is as 
functions of m^/fib- The superconductive state is gapless because the equation Cj(p) = 
has solutions for finite p; in other terms, differently from the CFL or the 2SC models, 
where the dispersion laws have the form IjlUI) , in the gCFL model there are gapless modes 
(this happens also in the similar model with two flavors, the gapless 2SC model |14|). 
This result can be understood as follows. The chemical potentials of the various quarks 
are in general different, so that in general one expects ^ 0. It follows that the generic 
quasiparticle dispersion law has the form 



(18) e(p) = ±6fi + y/(p - Jif + A 2 , 

which shows that, for appropriate values of Sfj, there are gapless modes. 

Before being accepted, the tentative ground state must pass a stability test. For gluons 
this means that the Meissner masses, obtained by the eigenvalues of the polarization 
tensor, must be real. In the CFL phase, with m s = the Meissner masses are degenerate 
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with the value 



(19) 




( 



3G 



11 




) 



> 0. 



However, for a non zero strange quark in the gCFL phase the masses are in general 
different, due to symmetry breaking, and some of them are actually imaginary 0. This 
means that the gCFL state cannot be the true vacuum state, even if its free energy is 
smaller than that of the normal state. 

4. Color Superconductivity and compact stars 

Color superconductivity might be realized in compact stars. This follows from the 
following considerations. The BCS critical temperature is given by 



and in QCD Abcs is expected to range between 20 to 100 MeV, see eq. GTty . This result 
must be evaluated in the context of the thermal history of pulsars. These compact stars 
are formed after a supernova explosion. The temperature at the interior of the supernova 
is about 10 11 K, corresponding to 10 MeV. Then the star cools very rapidly by neutrino 
emission with the temperature going down to ~ 10 9 -10 10 K in about one day. Neutrino 
emission is then the dominating cooling process for ~ 10 3 years. When the star reaches 
the temperature of ~ 10 6 K, it cools down due to X-ray and photon emission, so that in 
a few million years it reaches a surface temperature ~ 10 5 K. Therefore for the largest 
part of its existence a neutron star has T < T c , which implies the possibility of color 
superconductive states if fi is large enough. Since the star temperature is much smaller 
than the typical BCS energy gap (T n . s J Abcs ~ 10~ 6 — 10~ 7 ) one can assume that the 
compact star is effectively at zero temperature. 

We have shown that QCD favors the formation of BCS condensates in idealized cases, 
e.g. two or three massless flavors of quarks. However in realistic cases the three quarks 
have different masses and, as a consequence, different Fermi momenta. Let us estimate 
the order of magnitude of the scales involved in the description of a neutron star with a 
quark core. We begin with the simple case of a free gas of three flavor quarks, assuming 
that the up and down quarks are massless and the strange quark has mass m s |16) . We 
also assume that the weak interactions are in equilibrium. One gets for quark chemical 
potentials and Fermi momenta the results 



(20) 



T c = 0.57 A BCS 



fJ'U — A 4 — g^ei PF — Mil 7 

1 d 

1 . r- 
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where jj, is average chemical potential 

(22) fi = ^ (fi u + ii d + /j, s ) 

and n e the chemical potential of the electrons. Notice that 

(23) J2 ViNi + n e N e = fiN q - n e Q , 

i—u,d.s 

where 

(24) N q = £ Ni, Q = ^N u -±(N d + N s )-N e . 

i—u,d,s 

This result reflects the fact that the u quark has electric charge +2/3 e, the d and s 
charge — 1/3 e, and the electron — e. The chemical potential for the electrons is fixed by 
requiring electrical neutrality. This corresponds to the following condition for the grand 
potential O at zero temperature 

BD, 

(25) Q = 7F = °- 

For each fermionic species, omitting the volume factor, Q is given by 

(26) = \ r p 2 (E(p)-n)d P . 

n Jo 

In our case we get 

(27) n = — V / p\E i (jp)- lH )dp+— p 2 (p-fi e )dp, 
where 



(28) E u4 (p) = p, E s (p) = vV+Tnf . 

An analytical expression can be obtained by performing an expansion up to the order 
rrig/n 4 . One gets 

(29) , e « g 
and 



(30) 
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The baryon density is obtained as 

(31) ^ = -3^ = 3^ E 

i— u,d,s 

With the same approximation as before one finds 



(32) p B ~ ''" 



n 2 



1 / N 2 

1 / m 



1 

2 V M 



We note that densities in the core are of the order of 10 15 g/cm 3 , corresponding to a 
chemical potential of the order of 400 MeV. 

Let us discuss the range of values p, ~ 400 MeV of the average chemical potential, 
with a strange mass of the order 200-300 MeV (this is the effective density dependent 
strange quark mass. With m s = 300 MeV one finds /x e = 56 MeV with Fermi momenta 

(33) p u F = 365 MeV, p d F = 418 MeV, p" F = 290 MeV , 

and a baryon density about 4.4 times the nuclear matter density. With m s = 200 MeV 
the result is [i e = 25 MeV and 

(34) p F = 384 MeV, p% = 408 MeV, p s F = 357 MeV , 

and a baryon density about 5.1 times the nuclear matter density. To go to baryon 
densities relevant to the central core of the star, i.e. densities from 6 to 8 times the 
nuclear matter density, one needs to go to higher values of fi and lower values of m s 
where the difference among of the Fermi momenta is lower. This can be seen using the 
approximate expression for p e : 

„ m 2 s d ml 5m 2 

(35) Pf~V~ , 

with 



(36) P F -P F *P U F ~P F 



9 

mi 



Ap 



In conclusion the conditions for the formation of color condensates can be indeed realized 
in compact star, more exactly in their cores. The precise exact nature of the condensed 
phase is however still matter of debate. In fact, it is quite likely that the CFL is not 
favored, due to the non vanishing strange quark mass. Since the gCFL phase is unstable, 
one has to search for an alternative, A possibility is offered by the inhomogeneous color 
superconductivity that we now discuss. 
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5. Inhomogeneous color superconductivity: LOFF phase with two flavors 

The relevance of the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase [HI EI EDI EH 113 

follows from the fact that, for appropriate values of <5/i, pairing with non-vanishing total 
momentum: pi + P2 = 2q ^ is favored. It is a well known fact that if fermions 
of different chemical potentials interact in an attractive channel, pairing occurs only if 
the difference in chemical potentials 8^, does not exceed the limit Agcs/V^, called the 
Clogston-Chandrasekhar limit. For 8fi much larger than the Clogston-Chandrasekhar 
limit the favored state is the normal phase, but in a window of intermediate values of 
Sfi pairing with non-vanishing total momentum is favored, as first shown in the con- 
text of condensed matter in The LOFF phase is characterized by gap parameters 
not uniform in space that, in the simplest case, have the form of a single plane wave: 
A(r) = Aexp(2iq • r). Both the case of two and three flavours have been discussed in 
the literature. Let us begin with the case of 2 flavors , i.e. let us assume that the strange 
quark is decoupled due to its (rather) large mass. In nature flavor symmetry is broken 
not only explicitly by quark mass terms, but also by weak interactions. Therefore in the 
applications of the color superconductivity one has to take into account this symmetry 
breaking; for example in compact stars, considering only two flavors, isospin is broken 
by Sfi = fx u — /id 7^ 0, due to the process: 

(37) d —>■ u e" i> e . 

At the equilibrium (including the reaction u e~ — > d v e ) one has <5/i = — /i e . 

To simplify the problem I will consider the case of two massless quarks with chemical 
potentials \i u and fj, u given by 

(38) /i„ = fi + Sfi, \i d = /i - S/j, , 
where up and down refer to flavor. The condensate has the form 

(39) tyftf) ex ^\ l3 . 

Let Abcs be the value of the homogeneous BCS condensate. For Sfi < Abcs/V%, the 
energetically favored state is the homogeneous one. To show that for S/u, > Abcs/V% the 
Cooper pair might prefer to have non zero total momentum, we consider a four-fermion 
interaction modelled on one-gluon exchange, that is 

(40) L = -^7^A°^^A ^, 

8 

where A a are Gell-Mann matrices. In the mean field approximation it reduces to 

(41) L = -ie^tyfVfAe**' + c.c.) + (L^R), 
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where we have defined 

(42) r s e 2 ^ r = -i £ ^%(C^), 
and 

(43) A = gT s ■ 
The gap equation has the form 

. .gp . [dv f +s ,dE 1 

(44) A = t^-A / — / d£ 



4tt J_s s 2tt (E - p, + it sigmE) 2 - C - A 2 ' 

where /? = 4p 2 /ir 2 is the density of states, 5 the ultraviolet cut-off and £ the component of 
the quark momentum parallel to the Fermi velocity v, measured from the Fermi surface. 
Moreover 

(45) p, = 5[i — v • q. 
Performing the integration over the energy we get 

(46) l.ffi /S/ 5 9 ( £ -|»l), 



2 i 4^ y ^2 + A 2 

where 

(47) 0(e - \fi\) = 1 - 0(-e - fi) - 0(-e + fl) . 

It can be shown that there is a first order transition in Sp, between the homogenous 
state and the normal state and a second order transition between the LOFF state and 
the normal one. The first order transition occurs near 5p\ ~ A bcs I v2, while the second 
order phase transition is at the critical point 8p 2 = 0.754Ascs (see [8]). Near 5p 2 one 
has 



(48) A LOFF = ^l.7b7 5p 2 (8p 2 -5p) = 1.15 A BCS \ ^ ^ 

V ^BCS 

which is non vanishing, and the grand potential is given by 

(49) QLOFF - ^normal = -0.439 p(8p - 6 p 2 ) 2 , 

showing that the LOFF phase is indeed favored. 

In the more general case one can assume the ansatz 



(50) 



p 

A(r) = A e2iqm ' r 

rn — 1 
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(the previous case corresponds to P = 1). Only approximate solutions can be given 
for the gap 1)50(1 either near the second order phase transition, by the Ginzburg-Landau 
method, or by an effective gap equation that captures the main aspects of the structure 
represented by eq. (JSDJ. 

5*1. Gap equation in the Ginzburg-Landau approximation. - In the Ginzburg-Landau 
approximation the gap equation for condensate l|5U|) is as follows 

/ , 9Q 
< 51 > 0A-° 

with 

(52) ^ =P « A2 + ^ A 4 + 2 A6 
y ' p 2 4 6 

and 

a = — (1- 11(g)), n((7) = n(q,q), 

39 

2 P 

0= ^( C tk,q^,qm,q n )'5q k -q £ +q m -q„ , 

"P k,e,m,n=l 

2 P 

( 53 ) 7= ^2 / ^(qk,q^qm,qj,qi,qn)<$q k -q £ +q m -q j +q i -q„ • 

k,£,m,j,i^n—l 

Here <5(qk — q n ) means the Kronecker delta: S n> k and the following definitions have been 
adopted: 

(54) n (qi ,q 2 ) = +*f J d ± j^dt f^f[ME,5,,M), 



/j p-\-8 r~\-OQ j T7i 4 

(56) ^(qi,q2,q 3 ,q4,q 5 ,q6)= + J ^ d£ T^JlME, fa {<*})■ 
We have put w = vf w and 

(57) /;(£, {q}) = : ; 

E + it signS - 6» + (-l)'fc " 2 ELi • q k ] 
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moreover the condition J2k=i(~ l) fe< 3k = holds, with M = 2,4,6 respectively for II, 
J and K. The critical value 5fi2 is obtained by the condition that, at <5/j = <5/i2, a 
vanishes. This approximation is adequate to deal with most of the structures that can 
be obtained summing several plane waves. This analysis has been performed by Bowers 
and Rajagopal in [jj]. The most interesting case is offered by the cubic structures. They 
are formed either by six plane waves pointing to the six faces of a cube (the so-called 
body-centered-cube, bcc) or by eight plane waves pointing to the vertices of a cube (face- 
centered-cube, fee). In the former case the grand potential turns out to be bounded from 
below and smaller than all the other crystalline structures. In the latter case the grand 
potential is unbounded from below, because the coefficient 7 is negative. This represents 
a problem, because it shows that the analysis is incomplete and higher order terms, e.g. 
~ A 8 , should be included. 

5"2. Effective gap equation. - If one is far away from the second order phase transition, 
the Ginzburg Landau approach is not valid. One can use in this case an effective gap 
equation obtained by an appropriate average of the original lagrangian. The smoothing 
procedure is described in detail in jll| and we can recall here only the main points. The 
method employed uses effective quark fields where the large part of the quark momentum 
/rv has been extracted by a factor exp(i/iv) (v the Fermi velocity). Therefore in the 
condensate term of the lagrangian one has a factor 

(58) ^ ex P *(^« v « + MdVd + 2q m ) • r . 

m 

Next we multiply the lagrangian by some appropriate function g(r) and make an average 
over a crystal cell. In the gap equation the relevant integration momenta are small. 
Therefore one can assume that the fields are almost constant in the averaging procedure. 
Averaging the factor (|58|l produces as a result that the two velocities are antiparallel up 
to terms of the order of Sfi/n; moreover the term (I58|l is a substituted by its average 

p 

(59) A E (v,e )= £ A e //(v-n m ,4) , 

m— 1 

where £q is the quasiparticle energy and n m = q m /g. It can be shown that this procedure 
is valid for A not too small, i.e. far away from the second order phase transition. The 
exact form of the function A e ff depends on the function g(r) chosen for the average. In 
particular g(r) can be chosen in such a way that 

r A for (£,v) € PR 

(60) A e// =A8(E u )9{E d ) = I 

I elsewhere . 

where for each plane wave of wave number n m one has 



(61) 



E u4 = ±5(1 T gn m • v + VC 2 + A 2 
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and the pairing region (PR) is defined by the condition E u > 0, Ed > 0. This procedure 
simplifies the gap equation that assumes the form 

(62) PA = i^[ d ^[ d ^. 



2 J At: J 2n i% - l\ - A|(v,4) ' 

The energy integration is performed by the residue theorem and the phase space is divided 
into different regions according to the pole positions, defined by 



(63) e=^ + A|(v,e). 

Therefore we get 



(64) = §/X fe 4^^ ^ 2 + fc 2 A 2 
where the regions Pu arc defined as follows 

(65) P fc = {(v,£)|A B (v,e) = fcA} 
and we have made use of the equation 

(66) — =ln- 26 



gp Abcs 

relating the BCS gap Abcs to the four fermion coupling g and the density of states. 
The first term in the sum, corresponding to the region Pi, has P equal contributions 
with a dispersion rule equal to the Fulde-Ferrel P = 1 case. This can be interpreted as 
a contribution from P non interacting plane waves. In the other regions the different 
plane waves have an overlap. The free energy f2 is obtained by integrating in A the gap 
equation. At fixed S[i, fHs a function of A and q, therefore the energetically favored 
state satisfies the conditions 

on 

^ M ' -dq--^ 

and must be the absolute minimum. The result of this anlysis is that the body-centered- 
cube (bec) is the favored structure up to Sp, ~ .95Abcs- For larger values of <5/i < 
1.32 Ascs the favored structure is the face-centered-cube (fee). 
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6. LOFF Phase of QCD with three flavors in the Ginzburg Landau approx- 
imation 

We shall describe in this section some results for the the LOFF phase of QCD with 
three flavours, obtained in the Ginzburg Landau (GL) approximation We shall limit 
our presentation to the case of space modulation given by a single plane wave. The free 
energy per unit volume f2 in the GL limit is 

(68) si = a n + j2 [y^ + T A ' + £^T A ? A ') + °( A6 ) 

where fi n refers to the normal phase 

(69) n n = --^« + M 4 rf + ^) 



I2ir z I2ir z 
and one has assumed for the condensate the pairing ansatz 

3 

(70) < iPic C l5 i> pj >cx A/(r) e afn e ij7 

i=i 

with 

(71) Aj(r) = A J exp(2iq I T) . 

In l|69|) fij are chemical potentials for quarks or for the electron while the coefficients of 
(|68|l can be found in |12j . One works in the approximation of vanishing color chemical 
potentials fi^ = /Ltg = 0, which is valid for the present case since they vanish in the 
normal phase and we work near the secon order phase transition. (3— equilibrium is 
imposed together with the electric neutrality condition 

/ \ 9n 

(72) - = . 

djJL e 

These conditions, together with the gap equations, give, for each value of the strange 
quark mass, the electron chemical potential fi e and the gap parameters A/. Moreover 
one determines qi by solving, together with the gap equation and Eq. (|72l . also: 

(73) = I^A^+A^A^, 7=1,2,3. 

oqi dqi j-^ dqi 

The condition l|72|l gives 



(74) 
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which is valid up to terms of the order (l//x). This result is identical to the free fermion 
case, which was expected since one works near the transition point between the LOFF 
and the normal phase. It follows that 

(75) Sudu = S^us = 
and 

(76) dfids = 25/1 . 

To evaluate (|73|l ■ it is sufficient to work at the 0(A 2 ), which leads to q — 1.1997|5/Lt|. 

As to the orientation of qj , the results obtained in Ref. ^2] indicate that the favored 
solution has Ai = and therefore qi = 0. Furthermore q2 = q3 and A2 = A3. These 
results are consequences of the GL limit. In fact, as shown by Eqns. (|75() and l|76[) . the 
surface separation of d and s quarks is larger, which implies that Ai pairing is disfavored. 
On the other hand the surface separations of d-u and u-s quarks are equal, which implies 
that A2 and A3 must be almost equal. Finally q2, q3 are parallel because the pairing 
region on the u— quark surface at the GL point is formed by two distinct rings (in the 
northern and southern emisphere respectively), while for antiparallel q2, q3 the two rings 
overlap, which reduces the phase space available for pairing. 

The results one obtains can be summarized as follows. Leaving the strange quark mass 
as a parameter (though a complete calculation would need its determination, considering 
also the possibility of condensation in the qq channel) the free energy can be evaluated 
as a function of the variable m s / fi. One finds an enlargement of the range where color 
superconductivity is allowed, with the LOFF state having free energy lower than the 
normal and the gCFL phases. For example, for /i = 500 MeV, and with the coupling 
fixed by the value Ao = 25 MeV of the homogeneous gap for two flavors, at about 
m 2 / ' \x = 150 MeV the LOFF phase has a free energy lower than the normal one. This 
point corresponds to a second order transition. Then the LOFF state is energetically 
favored till the point where it meets the gCFL line at about m 2 /^, = 128 MeV. This is 
a first order transition since the gaps are different in the two phases. 

7. Stability of the LOFF Phase of QCD with three flavors 

In J7| the gluon Meissner masses in the three flavor LOFF phase of QCD where 
computed using the an expansion in l//i and the GL approximation. At the order of 
g 2 there is a contribution g 2 fi 2 /(2ir 2 ) identical for all the eight gluons; this result for 
the LOFF phase is identical to those of the normal or the CFL case as this term is 
independent of the gap parameters, see Q19f) . Another contribution to the polarization 
tensor is given by: 



(77) 



-Tr[iS(x,y)iH^iS(y,x)iH v b ] 
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where the trace is over all the internal indexes; S(x,y) is the quark propagator, and H£ 
is the vertex. The quark propagator S has components S** (i, j — 1,2), where each 
is a 9 x 9 matrix in the color-flavor space. At the fourth order in A one has: 



A, B = 1, • • • 9 are indexes in the basis A = (u ri d g , Sb, d ri u g , s r ,Ub, s g , db), A and A* 
are 9x9 matrices containing the gap parameters A/, pq is the energy, £ = |p| — /i, 
fi-A — P'A — M- S" 12 an d S 22 are obtained by the changes 11 <-» 22 and A <-> A*. 

In the region where the use of the GL expansion is justified, i.e. 128 MeV < m 2 /fj, < 
150 MeV all the squared gluon Meissner masses were found positive and therefore the 
LOFF phase of three flavor QCD is free from the chromo-magnetic instability. 

8. Neutrino emission by pulsars and the LOFF state 

As noted above neutrino emission due to direct Urea processes, when kinematically 
allowed, is the most important cooling mechanism for a young neutron star. When 
the temperature of the compact star is of the order of ~ 10 11 K neutrinos are able to 
escape, which produces a decrease in temperature. For smaller temperatures, e.g. below 
10 9 K, the direct nuclear Urea processes n — > p + e + v e and e~ + p — > n + v e , would 
produce rapid cooling. However they are not kinematically allowed, because energy and 
momentum cannot be simultaneously conserved. Therefore only modified Urea processes, 
involving another spectator particle, can take place. The resulting cooling is less rapid 
because neutrino emission rates turn out to be e v ~ T 8 , much smaller than the emission 
rate e v ~ T 6 due to direct Urea processes. 

These considerations apply to stars containing only nuclear matter. If hadronic den- 
sities in the core of neutron stars are sufficiently large, as we have seen, the central 
region of the star may comprise deconfined quark matter in a color condensed state. 
Therefore direct Urea processes involving quarks, i.e. the processes d — * u + e~ + v e 
and u + e~ —td + v e , may take place and contribute to the cooling of the star. In the 
color-flavor locked phase, in which light quarks of any color form Cooper pairs with zero 
total momentum, and all fermionic excitations are gapped, the corresponding neutrino 
emissivity and the specific heat C are suppressed by a factor e _A / T and the cooling 
is distinctly less rapid compared to quarks in the normal phase. However at densities 
relevant for compact stars the quark number chemical potential fi should be of the or- 
der of 500 MeV and effects due to the strange quark mass m s are relevant. In ^2] a 
calculation of neutrino emission rate and the specific heat of quark matter in the LOFF 



(78) 
(79) 



S 11 = So 11 + S n A [S 22 A* (S n + S n A S 22 A* S n )] , 
S 21 = S 22 A* (S n + So 1 A S 22 A* S n ) , 



where St? is the 18 x 18 matrix: 
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superconductive phase with three flavors was performed. It was shown that, due to the 
existence of gapless modes in the LOFF phase, a neutron star with a quark LOFF core 
cools faster than a star made by nuclear matter only. This follows from the fact that in 
the LOFF phase neutrino emissivity and quark specific heat are parametrically similar 
to the case of unpaired quark matter [e v ~ T 6 and C ~ T respectively). Therefore the 
cooling is similar to that of a star comprising unpaired quark matter. These results are 
still preliminary since the simple ansatz of a single plane wave should be substituted by 
a more complex behavior and it remains to be seen if a more complete calculation would 
confirm the results of [T8| . 

* * * 

I would like to thank R. Anglani, R. Casalbuoni, M. Ciminale, R. Gatto, N. Ippolito, 
M. Mannarelli and M. Ruggieri for a most fruitful collaboration on the themes covered 
by this review. 
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